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INVERSIONS RELATING STIRLING, TANH, LAH NUMBERS 
AND AN APPLICATION TO MATHEMATICAL STATISTICS 

GIACOMO DELLA RICCIA 

Abstract. Inversion formulas have been found, converting between Stirling, 
tanh and Lah numbers. Tanh and Lah polynomials, analogous to the Stirling 
polynomials, have been defined and their basic properties established. New 
identities for Stirling and tangent numbers and polynomials have been derived 
from the general inverse relations. In the second part of the paper, it has 
been shown that if shifted-gamma probability densities and negative binomial 
distributions are matched by equating their first three semi-invariants (cumu- 
lants), then the cumulants of the two distributions are related by a pair of 
reciprocal linear combinations equivalent to the inversion formulas established 
in the first part. 



1. Introduction 
The usual form of an inversion formula is 

i i 

where {a(n, m), A(n, m)} is a pair of inverse numbers and {f{n)}, {g(n)} are nu- 
merical sequences. An extensive study of inverse relations and related topics can 
be found in [xj . In this paper we will consider inversion formulas of a more general 
form: 

g(n, m) = ^ a(n, i)f(i, m) <-> f(n, m) = A(n, i)g(i, m), 

i i 
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with double-sequences {/(n, m)}, {g(n,m)} and corresponding number inverses 
{F(n, to)}, {G(n, to)}; in other words, we are interested in inverse relations convert- 
ing between arrays. Of course, from the general formula, one can get inverse rela- 
tions for sequences by fixing one of the two indices. As we shall see, the case m = 1 
is of particular interest because, in general, it is associated with identities involving 
important classical numbers. We used this approach to find relations connecting 
Stirling { [ ™J , {arctan t(n, to), tangent T(n,m)} and Lah L(n,m) numbers. 

Then, in analogy to the Stirling polynomials o~k(x) case, we introduced tangent 
5k(x) and Lah \k{x) polynomials and derived connecting relations. Most of the 
identities we obtained seem to be original. Finally, we discussed the above results 
in the light of a problem in Statistical Mathematics dealing with semi-invariants 
(cumulants) of shifted- gamma probability densities g(d; a, b, c) — r(?9 + c; a, b) 
and negative binomial distributions nb(m;r, A). We showed that if g($; a,b,c) and 
nb(w; r, A) are matched by equating their first three cumulants, then the cumulants 
7(71) and r]{n) of the two distributions are related by reciprocal linear combinations 
equivalent to the array inversion formulas established previously. 



2. Relations between Stirling, tanh and Lah numbers 

We will use Stirling ["] , {^} arctan t(n,m) and tangent T(n,m) (PQ, see 
pp. 258-260), and Lah L(n,m) (jS], see pp. 43-44) numbers, but with scale factors 
and appropriate notations: 

9{n, to) = (-l) (n - m)/2 i(n, to), 6(n, to) = (-1 ) ( "- m)/2 T(n, to); 
n 

TO 

_ 77 1 

l(n,m) = (-l) n - m L(n,m) = (-l) n L(n,m) = — 

ml 

The orthogonal relations satisfied by {9(n, m), 0(n, to)} and | , { ™j} | are of 
the form '^2 i a(n,i)A(i,m) = '^2 i A(n,i)a(i 1 m) = [m = n] and it is easily verified 
by direct calculation that these relations are also valid for {7(n, to), L(n, to)} . 

Recursions and egfs are simply obtained by introducing scales in the ordinary 
relations. Since the egfs of 9(n,m), 0(n,m) egfs are of particular interest, we give 
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the following explicit derivation. The known egfs pQ are 

t m (v) = > t(n,m)— = — -arctan m v, T m (v) = T(n,m)— = — -tan m z;. 
' n\ ml n\ m! 

n n 

Since t m (v), T m (v), as functions of v, and to have opposite parity, t(n,m) = 
T{n, to) = when n — to is odd. Thus, with n — m is even, putting v = iu, 
v 2 = —1, arctan(iu) = largtanhu and tan(iu) = itanhu, we get 

(2.1) M«) = E n «"" ro *(»»n»)^ =E„(-l)^^,^)f = ^argtanh ro n 

= ^i(5 ln T ± s) = t^t (Ej %+i (2j+ + i)i ) ; ^j+i = = (2j)!; 

e m (u) m (n) = £ n * n - m T(ra,m)£ = ^ n (-l)^T(n,m)^ = ^ tanh™ u 

- ^! U 2j + 1 (2j+l)! J ' U 2i+1 - \- L ) J 1 2j+l - 2j+2 

Thus, scale factors (— l)("~ m )/ 2 change arctan t(n,m), tangent T(n,m) numbers 
in arctanh 6{n,m), tank 0(n,m) numbers. For brevity, the pairs {[^Jil™}}) 
{9(n, to), Q(n, to)} and {T(n, m), L(n, to)} will be called Stirling, tanh and La/i num- 
bers, respectively. Scaled numbers basic properties are listed in Table 1. 

Theorem 2.1. Numbers in each pair, Stirling, tanh and Lah, convert between 
numbers in the other two pairs. 

2u-l 

Let tanh it = j^j, v = ^—^ — , then from the egfs listed in Table 1 we get: 

E„0(n,m)^ = ^tanh™ U = ^ (^) = ^L&m)^ 

(2.2) = Ej L(i, to)1 (^i) 1 = E t Hi, m) E„ {!} £ 
= E„ (Eii(i.m)T^)^. 

These equations imply 0(n,TO,) = E™=m -^(*> TO ){Tl anci j by dualities and inver- 
sions (Table 1), 6(n,m) = EtJih rnj\*\ , - Z(n, *)*(*, «*) and {^} = 
Ei ^( rl ? *)0(*j TO )- Hence, Lah numbers convert between Stirling and tanh num- 
bers. In Table 2 we listed the identities that are derived by use of inversions and/or 
dualities given in Table 1 . Since in the third and fourth inversion formulas Stirling 
numbers convert between Lah and tanh numbers, and in the fifth and sixth tanh 
numbers convert between Lah and Stirling numbers, the proof is complete. 
The basic structure connecting tanh and Stirling numbers is the following. 



4 
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Recurrence relations 

L m J Lm — 1J irnl 1 V m J Lm- 1J ImJ ' 

m={l}=[-=o], Ki=m=[-=o]- 

#(n + 1, to) = #(n, m — 1) + ra(ra — l)9(n — 1, to), 
6(n + 1, m) = 6(n, m — 1) — m(m + l)G(n, m + 1); 
6»(0, m) = 6(0, m) = [m = 0]; 9(n,0) = 6(n,0) = [n = 0]; 0(1,0) = 6(1,0) = 0. 
l(n + 1, to) = (n + m)l(n, to) + 7(n, to — 1), 
L(n + 1, to) = — (n + m)L{n, m) + L(n, m — 1); 
l(n, 0) = L(n, 0) = 0; (-l) m Z(0, to) = 1(0, to) = + [to = 0]. 
Duality laws and Orthogonal relations 

JnJ = (_l)«-"»^-mJ. 0( n> m ) = 9 (_ m; _ n ). /(^ m ) = (_1)«-™I(-TO, -T»). 

Ej a(n, i)A{i, to) = Ei A(n, i)a(i, m) = [m = n]. 
Exponential generating functions 

E„ ffi = & (h Mi + 2-)) m , En U£ = ^ (^) m 

M«) = En Hn, m)^ = ± argtanh™ u = ± (§ In 
©m(w) =E»e(».n*)^ = ^tanh m U 

Table 1 . Basic properties of scaled Stirling, tanh and Lah numbers 

Corollary 2.2. Tanh numbers are finite sums of multiples of Stirling numbers, and 
inversely 

g(»,«-fc) = Eto''q)( m 7 1 )M, 

L- J = E-=o(-l)^) CVM" k)_^ 

&in+ k, n) = Eto(-i) fc -^! rr) rr 1 ) o> 
r - Eto(-i) fc ^! rf ) (" + r x )0(" +*.«+<)■ 

These relations are obtained from the first two inverse pairs in Table 2 with 
k = n — to, i replaced by n — i and the use of Lah numbers explicit expressions. 
As we know (jj], see p. 418), Stirling numbers L^J , { x+ } can be viewed as 
polynomial in x. Thus, Corollary (|2.2(l implies that 6(x, x — k), Q(x + k, x) can also 
be treated as polynomials; these have the following properties. 

Proposition 2.3. If k = 2j > 0, then 6(x,x — k), &(x + k,x) are polynomials in 
x having degree 4^ = 3j and leading coefficient gj^ji , (— l) 3 3 j \^ , respectively. 
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LmJ 



= (-2)- m [:L {£}=2—»{£} 

(n, m) = (-l)("- m )/ 2 i(n, m), 9(n, m) = (-l)("- m )/ 2 T(n, m) 
I(n, m) = (-l) n i(n, m), I(n, m) = (-l) m L(n, m) 



0(n,m) =EiL(»,m){7} 
fl(n,m)=E i i(n,0'EI 



C}=Ej(i.m)e(n,i) 



L(n,m)=Ei [?]©(*,"») ~ e(n,m) = Ei{?}i(*,m) 



I(n,m)=Ei«(n,i)U} 
I(n,m) = Eie(i,Hra 



{^} = E<e(n J i)T(i > m) 
["]=Ei^,"i)X(n,i) 



Table 2. Conversions between Stirling, tanh, and Lah numbers 



The proof is by induction on k applied to tanh numbers recurrence relations 
(Table 1) written with k = n — m notations. Details of the proof are omitted 
because they are the same as those used by Gessel and Stanley ([2], see p. 25 ) in 
their study on Stirling numbers structure. 

The first few cases are the following. 

9(x,x-2) = -Q(x,x-2) = sin® 
9(x, x - 4) = ^ - 2%) ; e(x,x-i) = - s % B ft 1 ) - ^ * Hi) ■ 

As pointed out in the Introduction, the general inverse relations in Table 2 yield 
interesting results in the case of m = 1, essentially because 

[3 = (-2) n " m (n - 1)!; {f} = 2 n " m . 
0(2n+l,l) = (2n)!; 6(2n + 1, 1) = {-l) n T 2n+1 ; 6{2n, 1) = 0(2n, 1) = 0. 

The first general pair of inverse relations gives Stirling numbers identities: 

0, even n " 

(2-3) ElLi(-ir^!{ "} = - 2"- 1 =^i!9(n,i). 

I 0(n, 1), odd n i=i 

The identity when n is even was given for the first time by Lengyel ([BJ, see p. 7), 
whereas the identity for n odd is new. The third inverse pair yields: 

n n 

(2.4) 2"" 1 = 5^e(n,i)»! «-» n! = 5^0(n,i)2 i- \ 

i=l i=l 
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showing that tanh numbers convert between factorials and powers of 2. From the 
fourth we get an inverse pair: 

n n 

8(n, 1) - ^(fMJHO'-^i - 1)! «- (-2)'^ 1 (n - 1)! = £i(n,i)0(i, 1), 
»=i j=i 

extending an inversion formula n! = E"=i -M 71 ; i)2* £! <-> 2 n n! = E™=i ^( n > 

given by Lah (0, see p. 207). Finally, the fifth and the sixth pairs disclose original 

identities involving two out of the Stirling, tanh and Lah numbers: 

e(M) = ££=i{7}(-i) < - 1 i! - (-i) B - 1 n! = Er=iTae(i,i) 

(-l)»- 1 n! = Er =i e^,l)]f - (-2)»- 1 (n-l)! = Er=i^, !)£(»,*)• 

3. Stirling, tanh and Lah polynomials 

Since 9{x, x-k) = £*L ( El and x-k) = k\Q (V) vanish for 
x = 0, 1, . . . , k , £tm/i polynomials 5k(x) and La/i polynomials Xk(x) can be defined 
by rules similar to <Jfc(a;) = / a£±i [3] used for Stirling polynomials. For 

clarity, we recall that x^- = a; (a: - 1) . . . , (a; — k) = (-l) fe+1 (fc - 

Definition 3.1. 

£ fe (x) = <^> ~ ^(fc-,) = (-l)W^; 
A fc (x) = ^i ~ Afc (fc- 2 ;) = -Sf^ = (-l) fc + 1 A fe (x). 
In the case of integers n, m, the above definitions assume the form 

77,! 77-! 

(3.1) 9(n,m) = -. -— <5 n _ m (n); 6(n,m) = - 7 L -TT<5„- m (-m). 

(m — 1)! (m — 1)! 

Corollary l|2.2|l with x instead of n and a factor a;2±i divided out, yields: 

Proposition 3.2. Tanh polynomials Sk(x) are finite sums of multiples of Stirling 
polynomials ak(x), and inversely 

Sk(x) = J2t=o ^~i l )°k-i(.x - i) ~ 40 - x) = Ei= C^-I 1 ) 5 ^ - 
^) = E < fc =o( IB 7 1 )*fc-<( !C - < ) ^ ^(^-^) = E-=ote 1 )^- :c )- 

Proposition 3.3. The generating functions of xSk(x) and x8k{k + x) are: 
Efc xSk(x)u k = (ucothu) 21 = (^2j 2 2j ■^jt'" 2 "'^ > Bernoulli numbers; 
E fc x4(fc + .x)« fe - (iargtanhu) 1 = ( JL l n * = (V :: ^^) ". 
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The egf of 6{n + k, n) 1|2.1|) and definitions \ yield 

> 0(re + fc,n)- --— = / 7 TTjOkin + k)-. — — arg tanh u, 

^ v ; (n + fc)! ^(n-1)! + n! 



xSk(k + x)u k — ( —arg tanh u ) 



A' 

where, on the basis of the "polynomial" argument, we replaced n with x. Using the 
egf of Q(n + k, n) and proceeding as above, we obtain J^k xSk(x)u k = (ucothu) x . 

Proposition 3.4. Tanh polynomials 5k(x) satisfy the recurrence relation 

(x + l)S k (x + l) = (x- k)8 k {x) + {x- l)5 k - 2 {x - 1); xS Q (x) = 1. 
5k(x), k = 2 j > 0, has degree k/2 — 1 and 8k(x) = 0, k = 2j + 1. 

The recurrence relation is obtained by dividing out a common factor x— in the 
recurrence relation of 9(n,m) (Table 1), written with n = x and m = x + 1 — k. 
(Compare with (x + l)u^{x + 1) = (x — k)ak(x) + xak-i(%), (0 Exercise 18, 
Chapter 6)). 5k(x), k — 2j > has degree k/2 — 1 follows from Proposition 1|2.3(1 
and 5k(x) = 0, k = 2j + l, because 6(x, x— k) = Q(x+k, x) = 0, fe = 2j + l. (Recall 
that CTfc(x), k > 0, has degree k — 1, 0). The first few cases are the following: 

4(l) = 2 fc |f + [fc = l]; fc4(0) = 2 fe (2 fe -2)|f, k > 0; 

*fc(-i) = -^jr; M2i + i) = 5sipr 

<y 2 (s) = |; *4(a:) = ^k(V) - 3^; <5 6 (ir) = 3^31(2) - 3^5 (V) + 3^5- 
A companion expression of (|2.3|l . derived from Proposition l|3.2|) . is 

E^=i( 8 7 1 )^-i(x-<) = | 0, ° ddfc ; v k (x) = J2 ( x "^Sk^ix-i). 

I <5fe(a;), even fc i=i ^ ' 

Ejii ( a: 7 1 )^ r2 J+ 1 " i ( a ' ~ *) = is a new Stirling polynomials identity. 
The properties of x\k (x) follow at ounce from x\k (x) = . 

Proposition 3.5. Polynomials x\k(x) — (?) have degree k. Their generating 
function and recurrence relation are 

T, k >o x\ k {x)u k = (1 + u) x , Y, k >o xX k(k + x)u k = (j-r^rp, 

(x + l)\ k (x + 1) = (*+*) = ( fc * x ) + (*) = x [Xk-i(x) + Xk(x)} , xXo(x) = 1. 
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4. An application to a problem in Mathematical Statistics 

We now show that the above results have an application in a problem of Statis- 
tical Mathematics dealing with semi-invariants (cumulants) of shifted- gamma den- 
sities g($] a, b, c) = Y{fi + c; a, b) and negative binomial distributions nb(zu; r, A). 

For the sake of completeness, we first recall some standard definitions. 



g{&, a, 6, c) = r(i? + c; a, b) = < 



^(tf + c^cxph^], #>-c, 
otherwise 



^A) = T ^g TT ( T i I ) r ( T A I ) ro , A>0 

From the moment cgfs M Sfl (£) = e~ ct /(l - 6t) a , M„ b (t) = 1 /[l - A(e* - l)] r and 
the cumulant cgfs lnM ss (rj), In M„(,(i): 

lnM 8fl (t) - E„> 7HS = (-c + a6)i + E„>iaf i", 
In M n6 (t) = En>oV(n)% = rE m >o £(e* - 1)™ 
= E„> [E m >oKm - l)!OA™] £ [with (e* - 1)™ = m! E„> 

we get the nth cumulants 7(71), n(n) of a,b,c), nb(w; r, A): 

-c + afe, n = 1, 



(4.1) 7 (n) 



(n - 1)! afe™, n > 1; 



(4.2) «(n) = E^=iK™-l)!2 m -"OA™, " > °- 

Definition 4.1. Two distributions g{d; a,b,c) and nb(zn; r, A) are said to be 
matched if their first three cumulants are equal. 

By equating the first three cumulants 

1st cumulant (mean) = fi = — c + ab = rX 

2nd cumulant (variance) = a 2 = ab 2 = rA(l + A) 

3rd cumulant = 7 (3) = 2! ab 3 = 77(3) = 2! rA(l + A) (1/2 + A) 

we get matching conditions 



„ _ A(l+A) I _ 1 /o 1 \ _ _ ab _ rX . 
^ <^ " — ' (1/2+A)^ ' 0— i/^TA, C — 1+2fe - 1+2A , 

r = T^TT4> A = 6-l/2. 



For convenience, we will use scaled cumulants r/(n) = —n{n) and 7(71) = —7(71). 
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Lemma 4.2. In a matched pair {g{"&; a,b,c),nb(zu; r, A)}, cumulants 7(1.) and 
77(77) are polynomials in A having degree n 

= Em=i 7K m)X m , n > 0, 77(71) = Em=l w)A m ; 

(4.4) 7(71, m) = 2 m (m - 1)! [?(n - l,m- 1) + 2m7(n - l,m)] , 

(4.5) r ? (n,m) = 2 m (m-l)!{^}. 

The lemma holds for 77(71) as we already know 1)4. 2[l . It is also true for 7(77) 
because if {a, 6, c} are replaced by the matching values 14.3( 1. then from l|4.1|) : 



7(77) 



2A, 77 = 1 

2"(t7 - 1)!A(1 + A)(l/2 + A)"- 2 = E^=i 7(n, m)A™, 77 > 1, 



' n-2 * 



n-2 



X" 



A(l + A)(l/2 + A)«- 2 = E™ =1 Vn-W3— ' Vn-m-l,^ 

and 7(77, 777) = 2" l (m — 1)! [7(77 — 1, m — 1) + 2ml(n — 1, m)l , 77 > 1. 

Theorem 4.3. In a matched pair {g($; a,b, c),nb(w] r, A)}, cumulants 7(77) a?7<i 
77(77) 0/ i/ie i77;o distributions are related by reciprocal linear combinations 



7(77 + 1) = ^2 e ( n > + 1) 



77(77 



1) = ^0(n,i)7(i + 1), 77 > 0, 



i=0 



i=0 



where 0(77,777) = {-l) {n - m)/2 t(n, 777) and 6(77,777) = (-l) ( "~ m)/2 T(77, m). 

Since cumulants are polynomials in A, a solution #(77,777) and 6(n,m) exists if 
like powers of A on both sides are equal: 

(4.6) 7/(77 + 1, 777 + 1) = [A m+1 ] ^2 6(77, + 1) = 6(77, + 1, m + 1), 

i i 

7(77 + 1, 777 + 1) = [A m+1 ] °( n > *M* + 1) = £ 0(n ' l)r > ( * + !. m + !)• 
From (|4.5|) and Stirling recurrence relation, and JO}, we get 



(4.7) 



7/(77 + 1, 777 + 1) = 2™+ 1 m!{ ™+\} = 2™+ l m\ {£} + 2(m + l)^} 



7(« + 1, 777 + 1) = 2 m+1 m! [Z(z, to) + 2(to + l)7(i, m + 1)] . 
Substituting 1|4.7(1 into (|4.6() . rearranging the terms and dividing by 2 m m!, we obtain 
a recursion in m which automatically terminates after 77 — 777 recursive steps 



O - Elm 9(n, *)!(*, m) - 2(777 + 1) [{ m « J - EL m+1 0(77, 7)7(7 , 777 + 1) 

= 2 2 (m + 2)(m + 1) [{ m '| 2 } - Et™+i 9(n, 7)7(7, m + 2)1 =,...,= 0, 
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Thus {,^} -Y,i®(n,i)T(i,m) = 0, and 6(n,m) = (-l)("" m ^ 2 T(n, m) from The- 
orem (|2.1|) . Since 7(71, + 1) = 52™_ ^( n ' O'K* + 1) IS the reciprocal of rj(n + 1) = 
53" =0 0(n, 1)7(1 + 1), #(n, m) and 0(n, m) are necessarily number inverses, therefore 
9(n, m) = (— l)("~ m '/ 2 £(n, m). This completes the proof. Conversely, 

»7(n+l, m+1) = 0(n, 1)7(1+1, m+1) «-> 7(71+1, m+1) = 0(n, 1)77(1+1, m+1) 

implies Theorem Q2.lt and the inversion relations in Table 2, hence, </ie problem 
in Theorem 1(4. 3|) 0/ converting between cumulants is equivalent to the problem in 
Theorem l|2.1|l of converting between Stirling, tanh and Lah numbers. 

Finally, let us show that the inverse pair 1)2.4(1 : nl = X)"=o ^( n ' *)2 I_1 2™ _1 = 
^™ =0 0(n, «)«!, corresponds to an interesting particular case of Theorem 14.3(1 . In 
fact, let A — > while rX — a remains constant, then, as we know, the negative 
binomial distribution nb(w; r, A) tends towards the Poisson distribution p(w; a) 
with cumulants all equal to a. The matched shifted-gamma density tends towards 
3(1?; 4a, i, a), since, according to the matching conditions 1(4. 3() . a — > 4a, 6 — > 
|, c — > a, and its cumulants are obtained from (|4.1|) . Using the cumulant limiting 
values, one verifies that the reciprocal relations in Theorem (I4.3|l reduce to l|2.4() . 

5. Main results and Conclusion 

From general inverse relations converting between Stirling, tanh and Lah num- 
bers, we obtained a certain number of new identities by fixing m = 1 in the double- 
sequences involved. The same approach was used to study connections between 
ak(x) Stirling, 5k(x) tanh and A(x) Lah polynomials. Finally, we showed that the 
cumulants of a shifted-gamma probability density and a negative binomial distribu- 
tion can be related by reciprocal linear combinations fTheorem I4.3|) which turned 
out to be an instance of the tanh numbers inversion formula, hence, that this prob- 
lem and our problem fTheorem 12. 1(1 on number arrays, are equivalent. 
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